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Optical identification of the companion to PSR J1911-5958A, the 
pulsar binary in the outskirts of NGC6752 
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Abstract. We report on the identification of the optical counterpart of the binary millisecond pulsar PSR J 191 1-5958A, located 
in the outskirts of the globular cluster NGC 6752. At the position of the pulsar we find an object with V = 22.08, B - V = 0.38, 
U — B = -0.49. The object is blue with respect to the cluster main sequence by 0.8 magnitudes in B - V. We argue that the 
object is the white dwarf companion of the pulsar. Comparison with white dwarf cooling models shows that this magnitude and 
colors are consistent with a low-mass white dwarf at the distance of NGC 6752. If associated with NGC 6752, the white dwarf 
is relatively young, S 2 Gyr, which sets constraints on the formation of the binary and its ejection from the core of the globular 
cluster. 
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1. Introduction 

Recently, 5 millisecond pulsars have been discovered 
(D'Amico et al. 2002 1 towards the nearby galactic globular 
cluster NGC 6752. Three of them are located inside the 6'.'7 
core radius (Lugger et al. [1995 1 while the other two are out- 
side the 1 '92 half-mass radius (Diorgovski ll993> . at 2'7 and 
6'4, respectively. The latter of these is a binary millisecond 
pulsar, PSR J1911-5958A (hereafter PSR A), and has a low- 
mass (^0.19M o , assuming a 1.4M neutron star) companion 
in a 20 hour, highly circular (e < 10~ 5 ) orbit (D'Amico et al. 
2002 1. The pulsar period and period derivative suggest that it 
is a "canonical" recycled millisecond pulsar (see Phinney & 
Kulkarni 1994 for a review), and hence that the companion is 
likely a white dwarf. 

The large separation of PSR A from the cluster center is 
puzzling. Colpi et al. (|2002i have recently investigated possi- 
ble scenarios, and found that both for the case of a primordial 
binary and for an exchange or scattering event with other clus- 
ter stars, it is very difficult to explain both the pulsar's current 
position and its close circular binary orbit. Instead, they suggest 
that the binary may have been scattered to its current position 
by a binary composed of two 3-100 M Q black holes. 

One might learn more about the system's origin (and verify 
cluster membership) if one can confirm that the companion is 
a white dwarf and measure its mass and age. Therefore, we 



searched for the optical counterpart in archival data. We report 
here on the results. 



2. Observations and analysis 

We searched the ESO and Hubble Space Telescope archives 
for images coincident with the pulsar position (a^ooo = 
19 h ll m 42 s :7562, d 2 ooo = -59°58'26'.'900; D'Amico et al. 
2002 1. We found two images, taken with the Wide Field Imager 
(WFI) at the ESO 2.2 m telescope on La Silla. These observa- 
tions, 4 minute B and V-band exposures, were taken during the 
night of May 13/14, 1999. The seeing was poor, ~ 1'.'5 in V. 
However, both images showed a faint object at the pulsar posi- 
tion (see below). This object was also present in two HST ob- 
servations with the Wide Field Planetary Camera 2 (WFPC2; 
Holtzman et al. U995> . Both observations, U5FI07 and U5FI03 
(GO-8256), were imaged in the same filters and had similar ex- 
posure times, 42 s in F555W (hereafter V 555 ), 220 s in F439W 
(B439), 660 s in F336W (C/ 336 ) and 1 800 s and 1 693.5 s in 
F255W (nUV255) for the first and second field, respectively. 
The position of the pulsar coincides with the WF3 chip for 
the U5FI07 dataset, while it is on the WF4 chip on the other 
dataset. 



2.1. Astrometry 



Send offprint requests to: C.G. Bassa, 
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The WFI detector has an array of 8 CCDs (2 rows of 4), each 
CCD having a field of view of 8' x 16', a total of 33' x 34'. 
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The position of the pulsar was coincident with chip 6 of the V- 
band image. We found that there was some distortion over the 
whole chip. To minimize its effect we only used the upper half 
of this chip for the astrometric calibration. Stars on this 8' X 8' 
sub-image were compared against entries in the USNO CCD 
Astrograph Catalog (UCAC; Zacharias et al. 2OO0J. In total 68 
UCAC stars coincided with this image and their centroids were 
measured. Of these 43 were not saturated and appeared stellar 
and unblended. One outlier, having a total residual of 0732 was 
rejected. The remaining stars were used to calculate an astro- 
metric solution, fitting for zero-point position, scale and posi- 
tion angle. This solution has root-mean-square (rms) residuals 
of (X'06 in both right ascension and declination. The statistical 
uncertainty in the astrometry is thus 07084 in each coordinate. 

We used the astrometrically calibrated WFI image to obtain 
astrometric solutions for the two HST/WFPC2 datasets. First, 
the WFPC2 pixel positions were corrected for geometric dis- 
tortion and placed on a master-frame, using the prescription of 
Anderson & King (2003 1. We matched stars on the WFI im- 
age and fitted for zero-point position, scale and position an- 
gle against the WFPC2 master frame positions. Outliers having 
residuals larger than three times the rms residual of the fit are 
removed and a new solution is computed. This process is iter- 
ated until convergence. On average the converged astrometric 
solution used some 200-300 stars with rms residuals of the or- 
der of 0703 in both right ascension and declination. The final 
uncertainty in the tie to the UCAC system is dominated by the 
step from the UCAC to the WFI sub-image, and is ~ 0710 in 
each coordinate. 



2.2. Photometry 

We started with the pipe-line calibrated HSTfWFPC2 images, 
and used the HSTphot 1.1 (Dolphin 2000a I package for fur- 
ther reduction and photometry of the images. We followed the 
recommended procedures to mask bad pixels, defects, cosmic 
ray hits and hot pixels. Next we used the main task hstphot to 
find stars, measure positions, and determine calibrated photom- 
etry. The latter uses the aperture corrections, charge-transfer 
efficiency corrections and zero-points of Dolphin ( 2000b I. 



2.3. The counterpart to PSR J191 1-5958A 

The UCAC catalog is on the ICRS at the 20 mas level (Assafin 
et al. 2003 ). Including this uncertainty in the uncertainty of the 
astrometric tie, we obtain 95% confidence radii for the WFI and 
HST/WFPC2 frames of 072 11 and 07235, respectively. Within 
these radii there is a single object, see Fig.^ The HST/WFPC2 
positions and magnitudes are tabulated in Tabled 

The photometric measurements of the object in the two 
datasets are consistent for all filters except nUV255- The object 
is 0.7 magnitudes brighter in the first dataset than in the sec- 
ond. However, we compared the magnitudes of 240 stars that 
overlapped between the two datasets and found that, on aver- 
age, stars in the first dataset were brighter by 0.34 magnitudes 
in nUV255, which removes the observed discrepancy. We have 
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Fig. 1. Finding charts for PSR A. The upper image is a 1' x 1' 
subsection of a V-band image obtained on May 14,1999 with 
the Wide Field Imager on the ESO 2.2 m telescope at La Silla. 
The 2" tick marks indicate the position of PSR A. The lower 
image is a 20" x 20" sub-image of the HST/WFPC2 V555 
(F555W) observation. The 07235 (95% confidence) error cir- 
cle is shown at the position of the pulsar. 



not, however, found an explanation for this large offset between 
the two datasets. 

The magnitudes in the HST flight system filters (^336, 
B439 and V555) were transformed to the Johnson-Cousins UBV 
(Vega) system using the transformations by Holtzman et al. 
(1995). The resulting average magnitudes for the object in 
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Dataset 


Date (UT) 


R.A. (J2000) 


Decl. (J2000) 


nUV 255 


u 336 


6439 


V555 


U5FI87 


March 2, 2000 @ 10:29 


19 h ll m 42!756 


-59°58'26'.'87 


20.74 ± 0.23 


21.49 ±0.15 


22.40 ±0.21 


22.10 ±0.17 


U5FI83 


March 3, 2000 @ 13:54 


19 h ll m 42757 


-59°58'26'.'83 


21.44 ±0.24 


21.78 ±0.14 


22.49 ±0.17 


22.09 ±0.14 




U-B B-V 



Fig. 2. Color-magnitude diagrams of the two HST/WFPC2 
fields. The magnitudes are transformed from the HST/WFPC2 
flight system to Johnson-Cousins U BV using the prescription 
by Holtzman et al. ( 1995 1. The average colors and magnitudes 
of the counterpart to PSR A are indicated with error bars. The 
errors include the uncertainty in the distance modulus and red- 
denning. Also shown are Helium-core white dwarf cooling 
tracks (solid lines) by Serenelli et al. (2002 1. The models have 
masses of 0.183, 0.197, 0.230, 0.244, 0.300 and 0.336 M , de- 
creasing in mass to the red, for a metallicity of Z = 0.001. The 
age of the white dwarfs along the cooling track is indicated. 

this system are U = 21.96 + 0.10, B = 22.46 + 0.14 and 
V = 22.08 ± 0.11. The object is blue by about 0.8 magnitudes 
with respect to the cluster main sequence in B - V at the same 
B-band magnitude, Fig. [5] This is also confirmed by an instru- 
mental color-magnitude diagram constructed from the WFI B 
and V-band images. 

Given the density of objects on the ffi>T/WFPC2 chips, 
there is a ~ 1% probability of a chance coincidence in the 95% 
confidence HST error circle. However, there are only a few 
stars that are blue with respect to the cluster main sequence, 
which greatly reduces the probability of a chance coincidence. 
The blue color, together with the positional coincidence, gives 
us confidence that we have detected the companion of PSR A. 

3. Ramifications 

The minimum mass for the companion is constrained by the 
pulsar mass function. For a pulsar mass of 1 .35 M Q (Thorsett 
& Chakrabarty 1999 1 this minimum mass is 0.185 M . The 



lower limit increases for heavier pulsars, roughly by 0.004 M 
for every 0.05 M step in the pulsar mass. Assuming a random 
probability distribution for the inclination of the binary, we find 
that there is 90% probability that the companion mass is less 
than 0.5 M . Given this range of masses and the fact that PSR A 
is a recycled millisecond pulsar, it is likely that the companion 
is a Helium-core white dwarf. 

To verify whether our observations are compatible with 
a Helium-core white dwarf at the distance of the cluster, we 
compare our magnitudes with the predictions from the white 
dwarf cooling tracks of Serenelli et al. (2002i. We use their 
tracks for Z = 0.001, as this metallicity provides the best 
match to the metalli city of NGC 6752 ([Fe/H] = -1.43 + 0.04; 
Gratton et al. 2003 1. We also assume a V-band distance mod- 
ulus (m - M) v = 13.24 + 0.08 and reddening E B . V = 0.040, 
as recently determined by Gratton et al. (2003 >, and use the 
relative extinction coefficients listed by Schlegel et al. J1998I . 

Figure |5] shows the cooling tracks for Helium-core white 
dwarfs with Z = 0.001 for masses in the range of 0.183 to 
0.336 M . It appears that the magnitude and colors of the com- 
panion to PSR A are compatible with the two lowest-mass 
tracks, 0.183 and 0.197 M . Note that the 0.183 M Q model is 
below the minimum mass inferred from the pulsar mass func- 
tion. 

We have fitted the observed absolute UBV magnitudes 
against the predictions from the Z — 0.001 Helium-core white 
dwarf models. A x 1 statistic was computed for each entry in 
the model from the difference between the observed and mod- 
elled absolute magnitudes. Table[5]shows, for each model with 
a given mass, the properties of the white dwarf at the x 1 mini- 
mum. Both Fig.|2]as Table[2]show that the lowest mass models, 
0.197 to 0.244 M , are preferred. At these masses the white 
dwarf is rather hot, with 7/ eff a 1 1 000 - 16 000 K. 

Given these high temperatures, the counterpart is relatively 
young, 552 Gyr (see both Fig.[2]and Table[2}. The precise value 
strongly depends on the mass, since white dwarfs with lower 
masses have relatively thick hydrogen envelopes, where resid- 
ual hydrogen shell burning keeps the white dwarf hot. From 
the values listed in Table|2] one sees a jump in the cooling age 
between the 0.197 and 0.230 M models. This reflects a di- 
chotomy in the thickness of the hydrogen layer, where, above 
a certain critical mass, the thickness has been reduced by shell 
flashes early in the evolution (Althaus et al. 2001 1. As the shell 
burning is through the CNO cycle, the critical mass depends on 
the metallicity: -0.18 M for solar metallicity, -0.22 M for 
Z = 0.001 and -0.26 M for Z = 0.0002 (Althaus et al. EOPTl 
Serenelli et al. l2"U0H . 

We should note that the dichotomy in cooling age was not 
found by Driebe et al. (1998 1, who found the flashes hardly 
affected the thickness of the hydrogen layer. As a result, more 
massive white dwarfs cool slower and are older. However, for 
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Table 2. Fitting results for Helium-core white dwarf models 
with a metallicity of Z — 0.001. The observed absolute UBV 
magnitudes were fitted against the modelled values. The dered- 
dened, observed values U - B, B - V and My are given in the 
second row. 



Mass 


1 til 


* C 


IT] - R) n 


(R - V)n 


My 


Y 2 
Xy 


t M oJ 


(KJS.) 


(uyr) 


-U. J^t 


U. Jt 


o.o^t 




0.183 


10.6 


2.30 


-0.31 


0.08 


8.92 


0.9 


0.197 


11.7 


1.36 


-0.37 


0.03 


8.98 


1.1 


0.230 


14.7 


0.37 


-0.59 


-0.06 


9.14 


3.4 


0.244 


15.6 


0.27 


-0.64 


-0.08 


9.18 


4.2 


0.300 


19.1 


0.07 


-0.80 


-0.13 


9.29 


8.1 


0.336 


20.3 


0.05 


-0.85 


-0.14 


9.33 


9.5 


0.380 


21.6 


0.05 


-0.89 


-0.16 


9.34 


10.6 


0.390 


22.1 


0.04 


-0.91 


-0.16 


9.34 


11.1 


0.422 


23.8 


0.03 


-0.95 


-0.18 


9.39 


12.9 


0.449 


24.7 


0.06 


-0.98 


-0.19 


9.40 


13.7 



the mass range of interest here, their results are similar. Their 
0.195 M model has a cooling age of 1.2 Gyr, comparable to 
the age of the Serenelli 0.197 M model. For higher masses the 
cooling age decreases. The 0.300 M Q model by Driebe et al. 
( 1998 1 has a cooling age of 0.2 Gyr, down to 25 Myr for the 
heaviest model, with a mass of 0.414 M . 

Finally, we note the similarity between the position in the 
color-magnitude diagram of the companion of PSR A and the 
companion of PSR J0024-7204W (Edmonds et al. 120021 in 
47 Tucanae. On the basis of optical variability with the orbital 
period (3.2 h) and eclipses of the pulsar radio emission, the 
companion is argued to be a heated main sequence star. For 
PSR A, however, we can exclude this possibility: the spin-down 
luminosity is too low to lead to heating to the observed temper- 
ature. 

4. Discussion and conclusions 

We have detected the optical companion to the binary millisec- 
ond pulsar PSR J1911-5958A, which is located 6.'4 from the 
center of NGC 6752. The companion is blue with respect to the 
cluster main sequence by 0.8 magnitudes in B - V and com- 
parison of its colors and magnitude with white dwarf models 
shows that it is consistent with a Helium-core white dwarf at 
the distance of NGC 6752. 

Irrespective of the cooling models we use, we find that the 
white dwarf is at most 2 Gyr old. This age is similar to the 
<: 0.7 Gyr the binary can be expected to stay in the outskirts if 
it is currently on a highly eccentric orbit in the cluster (Colpi et 
al. 2002 ), which suggests that the white dwarf formed during, 
or shortly after, an encounter that also ejected the binary from 
the core, in an exchange interaction involving a binary with 
another star or binary. In the scenarios in which the binary was 
formed in the periphery, or scattered by a binary black hole, the 
coincidence of the two time scales has to be due to chance. 

One would expect the characteristic age of the pulsar to 
be similar to the cooling age of the white dwarf, as the pul- 
sar starts to spin down after the cessation of mass transfer, 
while the white dwarf starts to cool. This seems not to be 



the case for PSR A, for the characteristic age of the pulsar is 
P/(2P) ~ 17 Gyr (D'Amico et al. l2002> . It may be instead that 
the assumption underlying the characteristic age is wrong, and 
that the period at which the pulsar started spinning is similar to 
the current one. 

Compared to other white dwarf companions to millisec- 
onds pulsars, the companion to PSR J1911-5958A is bright, 
V = 22, which opens the possibility to determine detailed phys- 
ical parameters (e.g., Van Kerkwijk et al. 1996 Callanan et al. 
1998). For instance, one could measure the white dwarf tem- 
perature and surface gravity through spectroscopy. By compar- 
ison with models, this would lead to a mass and radius of the 
white dwarf, as well as a much more precise cooling age. The 
radius would allow one to confirm the association of the binary 
with NGC 6752. If the association is confirmed, the more accu- 
rate distance to NGC 6752 provides additional constraints on 
the radius and thus the mass of the white dwarf. Combining the 
white dwarf mass with a radial-velocity orbit of the white dwarf 
(and thus a mass ratio), would give the mass of the pulsar. 
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